Sorting by Periodic Potential Energy Landscapes: Optical Fractionation 
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Viscously damped objects driven through a periodically modulated potential energy landscape 
can become kinetically locked in to commensurate directions through the landscape, and thus can 
be deflected away from the driving direction. We demonstrate that the threshold for an object to be- 
come kinetically locked in to an array can depend exponentially on its size. When implemented with 
an array of holographic optical tweezers, this provides the basis for a continuous and continuously 
optimized sorting technique for mesoscopic objects called "optical fractionation" . 



Many natural and technologically important processes 
involve classical transport of small objects through mod- 
ulated potential energy landscapes. While the generic be- 
havior of modulated transport is well understood in one 
dimension 1], fundamental questions remain for higher 
dimensions. Colloidal particles flowing through arrays of 
optical tweezers 0, Q provide a uniquely accessible ex- 
perimental archetype for this class of problems. Experi- 
ments on transport through square arrays have revealed a 
Devil's staircase hierarchy of kinetically locked-in states 
as a function of orientation Within each locked-in 
state, particles select commensurate paths through the 
array independent of the driving direction. The ability 
to selectively deflect one fraction out of a flowing mixture 
was predicted to be useful for sorting and purifying 
mesoscopic samples. This Letter describes a practical 
implementation of this process, which we term optical 
fractionation. Examining the kinematics of optical frac- 
tionation further reveals that the underlying lock- in tran- 
sition can be exponentially sensitive to size. 

Optical fractionation exploits a competition between 
optical gradient forces exerted by an array of discrete 
optical tweezers Q and an externally applied force, as 
shown in Fig.Q A driven particle's trajectory can be de- 
flected enough by an encounter with one optical tweezer 
to pass into the domain of the next, and so on down the 
line. Such a trajectory is said to be kinetically locked-in 
to the array. Optical fractionation's selectivity emerges 
because objects with different sizes, shapes or composi- 
tions can experience substantially different potential en- 
ergy landscapes in the same optical intensity distribution, 
and the periodicity of the lattice emphasizes these differ- 
ences. More strongly driven or weakly trapped objects 
escape from the array and flow away in the driving di- 
rection. The two resulting fractions can be collected in 
separate microfluidic channels downstream. 

To demonstrate optical fractionation in practice, we 
studied the transport of water-borne colloidal particles 
flowing past a linear array of optical tweezers projected 
with the holographic optical tweezer technique 0. The 
colloidal suspension was confined to a 4 mm x 0.7 mm x 
40 /im glass channel formed by bonding the edges of two 
glass cover slips. Capillary forces at the channel's inlet 




FIG. 1: Principle of optical fractionation, (a) Different types 
of particles are driven by external force Fo through an array of 
optical traps inclined at angle with respect to Fo . Strongly 
interacting particles (ai) become kinetically locked-in to the 
array and deflected, while the others (02) are not. (b) Tra- 
jectories for large (ai = 0.79 /im) and small (02 = 0.5 /mi) 
spheres calculated with Eq. © for experimental conditions 
described in the text. 



were used to create a flow of about 60 /im/s along the 
midplane. This flow carried a mixture of monodisperse 
silica spheres of radius a\ — 0.79 /im (Duke Scientific 
Lot No. 24169) and 0,2 — 0.5 /im (Duke Scientific Lot 
No. 19057), both of which can be tracked to within 30 nm 
in the plane at 1/60 sec intervals using digital video mi- 
croscopy 5] . The two sizes can be distinguished reliably, 
and typical trajectories appear in Fig. O 

Colloidal silica spheres are roughly twice as dense as 
water and settle into a monolayer just above the lower 
glass wall of the channel, with the smaller spheres floating 
about 1 /im higher. Given the Poisseuille flow profile in 
the channel, the larger spheres travel somewhat slower, 
with a mean speed of u\ = 13 ± 2 /im/s, compared with 
the smaller spheres' U2 = 17 ± 9 /im/s. The spheres 
experience a viscous drag, Fj = jjtij, characterized by a 
drag coefficient, 7^, that depends both on their radii, a^, 
and on proximity to surfaces Q . 

Twelve discrete optical tweezers were arranged in a line 
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FIG. 2: Optical fractionation of bidisperse silica spheres, (a) Representative trajectories for a± = 0.79 /im at 1/60 sec intervals, 
(b) Trajectories for a 2 — 0.50 /im obtained simultaneously, (c) Time-averaged areal density m(f) for a± relative to the mean, 
no- Data compiled from 30,000 trajectories over 4 hours, (d) Simultaneously acquired data for 0.50 /im radius spheres compiled 
from 45,000 trajectories. The color bar indicates m/no for both data sets, and the scale bar denotes 10 /im for all four panels. 



with center-to-center spacing b = 3.6 ± 0.1 jam oriented 
at = 12.0° ± 0.5° with respect to the flow. Each trap 
was powered by 1.7 ± 0.8 mW of laser light at 532 nm, 
which slightly exceeded the empirically determined lock- 
in threshold for the larger spheres, given and b. 

The trajectories in Fig. Efa) and (b) demonstrate that 
the larger spheres are indeed systematically deflected by 
the array of traps, while the smaller spheres are not. Con- 
sequently, the array creates a shadow in the distribution 
of large spheres into which the small spheres can flow. 
Because the purification of small spheres and concentra- 
tion of large results from lateral deflection, this optical 
fractionation process can proceed continuously, in con- 
trast to most competing techniques Q. 

These qualitative observations can be made more com- 
pelling by considering the statistics of tens of thousands 
of trajectories compiled in Fig. Efc) and (d). Here, 
we plot the two populations' time-averaged areal den- 
sities rij(r) normalized by their means. The separa- 
tion's quality is assessed in Fig. [31 with Q(f) = [ni(r) — 
n 2(f)} / [n\(f) ^^(f)}, which reaches unity in regions con- 
taining only large spheres, and minus one in regions with 
only small spheres. A transverse section taken along line 
A in Fig. UJa) and plotted as squares in Fig. Efb) reveals 
a thoroughly mixed sample with Q(h) = approaching 
the traps. A similar section along line B downstream of 
the array demonstrates roughly 40 percent purification of 
both large and small spheres. Much of the background 
can be attributed to the traps' nonuniform intensities Q, 
with large particles tending to escape from the weakest 
traps. In denser suspensions, this escape rate is increased 
by collisions. Both processes can be mitigated by project- 
ing multiple lines of traps |4| . 

Optical fractionation's ability to distinguish objects 



arises as a general and previously unappreciated fea- 
ture of transport through periodically structured envi- 
ronments. Analyzing such transport not only provides 
insights into optimizing practical sorting, therefore, but 
also sheds new light on a range of analogous processes. 

The potential energy landscape presented by an optical 
trap array is a convolution of the traps' intensity profile 
1(f) with an object's optical form factor f(f): V(f) = 
—1(f) o f(f). The total applied force then is 

F = V[/(r)o/(r)]+Fo, (1) 

where Fo is the driving force. Colloidal particles' mo- 
tions are overdamped, so that the resulting velocity is 
v = F/j. The associated trajectory has a component 
v x along the row of traps and another v y perpendicular. 
Although overdamped transport through a periodic po- 
tential is well understood for one-dimensional systems |l[ , 
few analytic results are available for the inclined line, and 
fewer still incorporate thermal or quenched randomness. 
Consequently, we focus on the kinematic limit in which 
both driving and trapping forces exceed random thermal 
forces so that trajectories may be treated deterministi- 
cally. We then estimate the threshold for an object to 
escape from an array of optical traps, and thereby estab- 
lish the selectivity of optical fractionation. 

Any trajectory entrained by the traps, such as the ex- 
ample in Fig. mb), is characterized by turning points 
where v y = 0. Conversely, any trajectory without such 
turning points must be unbounded. This establishes as 
the maximum possible locked-in deflection angle 

sm0 m < max{9, [-1(f) o f(f)]} /F . (2) 

At a given 0, objects are either deflected or not with 
a selectivity set by the dependence of sin m on material 
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FIG. 3: (a) Spatially resolved quality of separation Q(f) ob- 
tained with a single line of optical traps. The cross-section 
transverse to the flow direction along line A is plotted as 
squares in (b) and provides a baseline profile for the suspen- 
sion's composition before optical fractionation. The cross- 
section along B, plotted as circles in (b) shows the influence 
of the trap array. The curve in (b) is a guide to the eye. 



properties. To estimate # m , we model the array as a 
periodically modulated line of light with intensity Iq: 



1(f) = I A(y) ^2 a n cos(n/cx), 



(3) 



n=0 



where k = 27r/6, the dimensionless transverse distribu- 
tion A(y) is symmetrically peaked around A(0) = 1, 
and the coefficients a n account for the tweezers' detailed 
structure with a n = 1. 

Convolving first along the x direction by applying the 
Fourier convolution theorem to each term in the sum, 
and then noting that cos nkx < 1 yields 



< — max ' 
~ F 



A (y) o^2a n f(nka,y) 



n=0 



(4) 



where Vb oc Io is the potential wells' depth and / is 
the form factor's Fourier transform along the x direc- 
tion. The array's periodicity thus selects a discrete set 
of wavenumbers from the continuous / whose depen- 
dence on a endows optical fractionation with exceptional 
size selectivity. This is most clearly demonstrated if / 
can be factored into inline and transverse components, 
f(nka,y) = f x (nka)f y (y). In this case, 



sin ( 



< q(a) ^2a n f x (nka), 

n=0 



(5) 



with q(a) = k(o)Vq/Fq and k(cl) = 
max{— d y [A(y) o f y (y)]}. Equivalent results can be 
obtained when / is not separable, but require a lengthier 
term-by-term analysis. 

We turn our attention first to the transverse contri- 
bution. If a particle is comparable in size to the op- 
tical tweezers' width, wq, then n(a) depends no more 
strongly on size than 1/a. For example, if A and f y 
are Gaussia ns of wi dths wo and a, respectively, then 
ft (a) oc + a 2 . Similarly, the potential depth Vo 

and driving force Fo scale as simple powers of a in the ab- 
sence of nonlinear effects such as Mie resonances 8] . Con- 
sequently, the prefactor q(a) describes no more than alge- 
braic sensitivity to size and material properties. Compa- 
rable selectivity is offered by other techniques such as gel 
electrophoresis and field flow fractionation |g , and would 
be obtained with an unmodulated line of light (qq = !)• 

The in-line contribution is more interesting. Because 
a particle's form factor vanishes outside the interval 
x G [—a, a], its Fourier transform depends very strongly 
on wavenumber. For example, a uniform dielectric cube 
aligned with the array has a separable form factor, with 
fx{x) = 1 for x G [—a, a]. Its Fourier transform, 



f x (nka) 



sin nka 



nka ' 

is bounded by the leading-order cumulant expansion 

n 2 k 2 a 2 



f x (nka) = exp 



6 



(6) 



(7) 



for ka < 7r. The equivalent result for a sphere [8| 
with f(f) = a/1 — r 2 /a 2 for r G [—a, a] is f x (nka) w 
(7r/2)exp (-n 2 /c 2 a 2 /6), with f y (y) w exp(-n 2 k 2 y 2 /6). 
Any smooth, bounded, positive-definite f(f) on x G 
[—a, a] similarly would surpass exponential sensitivity for 
ka > 1. Substituting Eq. into Eq. (0) therefore estab- 
lishes optical fractionation's exponential size sensitivity 
for 1 < ka < 7r. 

Comparable sensitivity to control parameters is ob- 
served in analogous transitions between sub-harmonic 
steps in drive charge density waves |£j and between ki- 
netically locked-in states in two-dimensional optical trap 
arrays [4| . Similar results also can be obtained for arrays 
of potential barriers, suggesting that arrays of optical 
tweezers also should be effective for sorting absorbing and 
reflecting particles that are repelled by laser light. This 
analysis also carries over to filtration by arrays of micro- 
machined posts 0, which therefore should be able to 
resolve objects substantially smaller than the inter-post 
separation under the right circumstances. 

Both f(nka,y) and the coefficients a n fall off rapidly 
with index n. Consequently, the sum in Eq. © is domi- 
nated by the first term, n = 1. This improves the approx- 
imations used in deriving Eqs. Q and Q and suggests 
that the result may be considered an estimate for sin 9 m 
rather than simply a bound. 
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To demonstrate this, we apply this analysis to our 
present experimental data, modeling the individual opti- 
cal traps as Gaussian potential wells 

no = w|>p(-^f), (8) 

with a 2 (a) « Wq + a 2 In this model, 

sin(9 m < q(a) exp (-^2) > ( 9 ) 

where g(a) = (2/y/e)Vo/(aFo). The weakest trap's oc- 
cupancy, 77, j, is inversely proportional to the particles' 
minimum speed, v m i n , as they pass through. Conse- 
quently, we can estimate the relative trap strength from 
the data in Fig. as q(a) = 2(1 — v m [ n /u). Similarly, 
the separation between the depleted region ahead of the 
traps and the position of maximum occupancy is 2 a (a). 
From the histograms in Fig. |5fc) and (d), we obtain 
q(a) = 1.6 ±0.1 and 0.9 ± 0.2, and a (a) = 0.85 ± 0.07 /im 
and 0.58 ± 0.07 jam for the large and small spheres, re- 
spectively [l]|. These results suggest 9 m = 14° ±1° > 
for the large spheres and m = 3° ± 2° <C for the small, 
which is consistent with the observation that only the 
large spheres are systematically deflected. 

The threshold sin 9 m depends only linearly on Vo and 
Fq. Thus, imperfections in practical trap arrays and fluc- 
tuations in the driving force can be largely compensated 
for by the substantially stronger dependence on particle 
size. Indeed, Figs. and |3| demonstrate robust size sep- 
aration despite more than 20 percent variation in flow 
velocity over the course of the experiment. 

Equation also offers insights into applying optical 
fractionation to nanometer-scale objects. Stokes drag 
scales linearly with a and the optical trapping poten- 
tial for Rayleigh particles scales as a 3 , so that q(a) oc a 2 . 
Sorting proteins or nanoclusters, therefore, will require 
enhancing Vo by four orders of magnitude. This might 
be accomplished by increasing the light's intensity, reduc- 
ing its wavelength and exploiting resonances. Even 
then, only algebraic size sensitivity should be expected 
for objects with a«A because ka <C 1 in this limit. 

We have focused on effects due to deflection transverse 
to the optical axis. Multi-dimensional separations could 
take advantage of Bessel beams' ability to exert con- 
trolled axial forces [l3| to distribute objects both trans- 
verse to and along the optical axis. 

In summary, we have demonstrated optical fractiona- 
tion for a model system of bidisperse colloidal spheres. 
This approach lends itself to continuous, rather than 
batch-mode fractionation, with continuous tuning and 
dynamic optimization over the entire size range acces- 
sible to optical trapping, i.e. nanometers to microme- 
ters. The abrupt transition from free flow to kinetically 
locked-in transport should offer exponential size selectiv- 
ity for objects larger than roughly 100 nm. Separation on 



the basis of other characteristics also can be optimized, 
although exponential sensitivity should not be expected 
in general. Our analysis focuses on the kinematic limit, 
Fob > Vo > ksT, which is both tractable and appropriate 
for weakly-trapped micrometer-scale colloid. Stronger 
trapping would require a more detailed treatment of ther- 
mally assisted hopping [l|. A substantially more so- 
phisticated analysis also would be required for higher- 
dimensional arrays. Similarly abrupt transitions should 
occur with variation of driving or trapping strength for 
vortex creep through patterned type-II superconductors 
[l^ , electron transport through two-dimensional electron 
gases [l5|? an d electromigration on crystal surfaces, with 
potentially useful applications resulting in each case. 
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